Genetic variation and structure of endemic and endangered wild celery (Kelussia odoratissima Mozaff.) quantified using novel microsatellite markers developed by next-generation sequencing.

Kelussia odoratissima Mozaff. (Apiaceae) is a native plant that has been traditionally consumed in Iran's food and pharmaceutical industries. Overharvesting of the taxon, especially at the beginning of the growing season, due to its considerable medicinal and economic value, is believed to be the main reason for the extirpating of this plant. The consequences of the severe anthropogenic impacts on the genetic diversity of populations are poorly known. In order to investigate the level of genetic variation and patterns of the genetic structure of K. odoratissima, we developed novel microsatellite markers using the 454 Roche next-generation sequencing (NGS) platform for the first time. Out of 1,165 microsatellite markers bioinformatically confirmed, twenty-five were tested, of which 23 were used to screen genetic variation across 12 natural populations. Our results showed that the average number of alleles per locus and the polymorphic information content (PIC) were 10.87 (range 7 to 27), and 0.81 (range 0.67 to 0.94), respectively. The mean observed and expected heterozygosities (± SD) across all populations were 0.80 ± 0.31 and 0.72 ± 0.14, respectively. The average pairwise FST among the populations was 0.37 (range 0.04 to 0.81). Bayesian and distance-based clustering, and principal coordinate analyses revealed at least four major genetic clusters. Although high level of structure can be explained by landscape topography and geographic distance, presence of admixed populations can be associated to seed or pollen dispersal. Contrary to expectations, the high level of genetic variation and lack of inbreeding suggest that overexploitation has not yet significantly purged the allelic variability within the natural populations in protected areas.

Does urbanisation lead to parallel demographic shifts across the world in a cosmopolitan plant?

Urbanisation is occurring globally, leading to dramatic environmental changes that are altering the ecology and evolution of species. In particular, the expansion of human infrastructure and the loss and fragmentation of natural habitats in cities is predicted to increase genetic drift and reduce gene flow by reducing the size and connectivity of populations. Alternatively, the 'urban facilitation model' suggests that some species will have greater gene flow into and within cities leading to higher diversity and lower differentiation in urban populations. These alternative hypotheses have not been contrasted across multiple cities. Here, we used the genomic data from the GLobal Urban Evolution project (GLUE), to study the effects of urbanisation on non-adaptive evolutionary processes of white clover (Trifolium repens) at a global scale. We found that white clover populations presented high genetic diversity and no evidence of reduced Ne linked to urbanisation. On the contrary, we found that urban populations were less likely to experience a recent decrease in effective population size than rural ones. In addition, we found little genetic structure among populations both globally and between urban and rural populations, which showed extensive gene flow between habitats. Interestingly, white clover displayed overall higher gene flow within urban areas than within rural habitats. Our study provides the largest comprehensive test of the demographic effects of urbanisation. Our results contrast with the common perception that heavily altered and fragmented urban environments will reduce the effective population size and genetic diversity of populations and contribute to their isolation.

Archaeological evidence of resource utilisation of the great whales over the past two millennia: A systematic review protocol.

Archaeological faunal remains provide key insights into human societies in the past, alongside information on previous resource utilisation and exploitation of wildlife populations. The great whales (Mysticete and sperm whales) were hunted unsustainably throughout the 16th - 20th centuries (herein defined as the modern period) leading to large population declines and variable recovery patterns among species. Humans have utilised whales as a resource through carcass scavenging for millennia; however, increasing local and regional ethnographic and archaeological evidence suggests that, prior to the modern period, hunting of the great whales was more common than previously thought; impacts of earlier hunting pressures on the population ecology of many whale species remains relatively unknown. Hunting guided by traditional ecological knowledge may have been sustainable and likely originated in societies that also incorporated opportunistic use of stranded individuals. The collation of georeferenced zooarchaeological data of the great whales between the 1st - 20th centuries CE worldwide will provide insight into the timescale and distribution of resource utilisation of the great whales and how this varied within and between societies, and may have changed over time. By comparing regions of known resource utilisation and breeding and feeding grounds of current-day whale populations, this information will subsequently be used to infer regions where whale populations were possibly lost or extirpated prior to detailed historical records. This systematic review protocol also provides a template for archaeologists, ecologists, and historians interested in using faunal remains to infer historical ecology and resource use of wild animal populations. The transparency of our data collection approach provides opportunities for reproducibility and comparability with future datasets.

Evolutionary and phylogenetic analyses of the barcoding region suggest geographical relationships among Blastocystis sp., ST3 in humans.

Blastocystis sp., has 21 distinct subtypes of which ST3 thought to be the most prevalent subtype. This study aims to analyze the global variations of ST3. In total, 496 sequences with more than 400 nucleotides from Asia, Europe, Africa, and America were included in this study. Results show that allele 34 was the most prevalent allele in all continents. The lowest and highest allele diversity were observed in Europe and Africa, respectively. The nucleotide diversity ranged from 0.0077 in Europe to 0.02 in Africa, and haplotype diversity ranged from 0.461 in America to 0.6 in Africa. The haplotype network and Bayesian structure showed at least two major clusters including Asia and Europe-Africa-America. Tajima's D values for all continents were negative and statistically significant, indicating an excess of rare nucleotide variants. Similarly, the Fu's FS test showed negative values for all regions, indicating an excess of rare haplotypes. Pairwise FST exhibited a high genetic differentiation between Asia and other continents. Mismatch analysis for all populations showed a unimodal distribution. Our findings indicate that there are two probable major clusters of Blastocystis sp. ST3, a cluster which is shared between Europe, Africa, and America, and a cluster which is restricted to Asia.

Taxonomy, Population Structure and Genetic Diversity of Iranian Leishmania Strains of Cutaneous and Visceral Leishmaniasis.

BACKGROUND: Despite the broad distribution of leishmaniasis in Iran, there is a little genetic information about the causative agents and epidemiological status of the disease. Genetic diversity of the parasite is suggested to be one of the factors, which influences the clinical manifestations of the disease. In this study, we investigated the genetic variations, population structure, and evolutionary history of Leishmania species from endemic foci of Iran. METHODS: Fifty-two isolates from humans, canines, and rodents from different endemic foci of Iran were used to sequence the N-acetyl glucosamine-1-phosphate transferase (Nagt) gene. Phylogenetic and structure analyses were performed to investigate inter- and intra-species diversity of the Leishmania isolates. RESULTS: In total, 10 haplotypes including L. major (n = 6), L. tropica (n = 2), L. infantum (n = 1) and L. turanica (n = 1) were identified across 52 isolates. Haplotype diversity (Hd) ranged from zero for L. infantum and L. turanica to 0.78 ± 0.136 for L. major. This study identified population structure of Leishmania isolates from different geographical regions of Iran. The results of the phylogenetic tree showed 4 distinct clades for each species of Leishmania. In addition, the highest intraspecies diversity was observed among L. major isolates. No correlation was observed between species and geographic distribution of haplotypes. CONCLUSIONS: Leishmania isolates were identified at the species level using the Nagt gene, low variation within species indicates conservation of this gene in Leishmania. The results provide knowledge into the evolutionary history of Iranian Leishmania isolates.

The taxonomic significance of ddRADseq based microsatellite markers in the closely related species of Heracleum (Apiaceae).

Many studies on Heracleum have shown poor correspondence between observed molecular clusters and established taxonomic classification amongst closely related species. This might reflect both unresolved taxonomy but perhaps also a lack of good genetic markers. This lack of appropriate and cost effective species-specific genetic markers hinders a resolved relationship for the species complex, and this in turn causes profound management challenges for a genus that contains both endemic species, with important ecological roles, and species with an invasive potential. Microsatellites are traditionally considered markers of choice for comprehensive, yet inexpensive, analyses of genetic variation, including examination of population structure, species identity, linkage map construction and cryptic speciation. In this study, we have used double digest restriction site associated DNA sequencing (ddRADseq) to develop microsatellite markers in Heracleum rechingeri. Genomic DNA from three individuals were digested with Sbf1 and Nde1 and size selected for library construction. The size-selected fragments were sequenced on an Ion Torrent sequencer and a total of 54 microsatellite sequences were bioinformatically confirmed. Twenty five loci were then tested for amplification, resulting in 19 of these being successfully amplified across eight species, comprising both the so-called thick-stemmed species (H. persicum, H. rechingeri, H. gorganicum and H. lasiopetalum), and thin-stemmed species (H. anisactis, H. pastinasifolium and H. transcaucasicum). Both Bayesian and distance-based clustering, and principal coordinate analyses clearly separated these into two groups. Surprisingly, three H. pastinacifolium populations were not separated from populations of the morphologically similar endemic species, H. anisactis, suggesting lack of genetic differentiation. Likewise, high genetic similarity was found between H. persicum and H. rechingeri populations, questioning taxonomic separation at the species level between these taxa. Further analyses are needed to re-evaluate the taxonomic significance of observed morphological variability currently applied to distinguish these sister taxa. Nevertheless, our results represent progress in the effort to develop cost-efficient molecular tools for species discrimination in this genus.

Comparative Analyses of Phytochemical Variation Within and Between Congeneric Species of Willow Herb, Epilobium hirsutum and E. parviflorum: Contribution of Environmental Factors.

The plants in the Epilobium genus are considered to have several important medicinal properties due to their unique chemical composition. Although metabolic profiles of medicinal plants are mainly controlled by genetic factors, their production is also to some degree influenced by environmental factors, thus, variations in the levels of phytochemicals may represent long-term ecological and evolutionary interactions. In order to depict the magnitude of natural variation in level of chemical compounds among conspecific populations of Epilobium hirsutum (n = 31) and E. parviflorum (n = 16), metabolite profiling of aerial parts of plants was performed with gas chromatography/mass spectrometry analysis. Putative identification and structure annotation revealed the presence of 74 compounds including 46 compounds considered secondary metabolites categorized into flavonoids (n = 8), phenolic acids (n = 26), steroids (n = 3), and terpenes (n = 5) across all populations. Although there was a considerable natural variation among conspecific populations, principal component analysis revealed a clear separation of populations of each species based on the second main principal component which was highly correlated with eight secondary metabolites. The level of secondary metabolites was significantly correlated between species (r = 0.91), suggesting shared metabolic pathways underlying the production of chemical compounds. In addition, redundancy and variance partitioning analyses by including bioclimatic variables and altitude revealed a significant contribution of elevation in explaining the total variation of secondary metabolites in E. hirsutum. Two-thirds of all secondary metabolites were significantly correlated with altitude in E. hirsutum. The large-scale geographic analyses of populations revealed additionally detected flavonoids and terpenes (E. hirsutum and E. parviflorum) and steroids (E. hirsutum) for the first time. This study provides significant information on additional chemical compounds found across the distribution range of the two ecologically important species of willow herb and emphasizes the importance of geographic-wide sampling as a valuable strategy to depict intraspecific and interspecific variability in chemical traits.

Large Scale Geographical Mapping of Essential Oil Volatiles in Heracleum (Apiaceae): Identification of Novel Compounds and Unraveling Cryptic Variation.

Conspecific populations of plants in their native range are expected to show considerable variation due to long-term ecological and evolutionary factors. We investigated the levels of secondary metabolites in Heracleum including H. persicum a valuable medicinal plant to depict the magnitude of cryptic variation and the potential significance of novel chemical traits. The essential oil volatiles from fruits of 34 populations from different species of Heracleum in Iranian distribution range and a native of H. sphondylium and an invasive population of H. persicum from Norway were analyzed with GC/MS. Out of 48 compounds identified, a contrasting pattern in the level of two major compounds, octyl acetate and hexyl butyrate was found among all studied species. Interestingly, a significant geographic pattern was observed; the hexyl butyrate/octyl acetate ratio was high (range 1.8 - 3.2) in the northwestern Iranian populations of H. persicum compared to that in northern and central populations (range 0.3 - 0.9). Four populations from Zagros mountains also exhibited a unique composition. Anethole was found in two populations of H. persicum from central Zagros, which has not been previously reported for essential oil of fruits of Heracleum so far. The results suggest high efficiency of large scale sampling from distribution range of species in identifying novel compounds. The unique pattern of geographic structuring also provides novel information to unravel cryptic variation in Heracleum.

Euglossine bees mediate only limited long-distance gene flow in a tropical vine.

Euglossine bees (Apidae: Euglossini) have long been hypothesized to act as long-distance pollinators of many low-density tropical plants. We tested this hypothesis by the analysis of gene flow and genetic structure within and among populations of the euglossine bee-pollinated vine Dalechampia scandens. Using microsatellite markers, we assessed historical gene flow by the quantification of regional-scale genetic structure and isolation by distance among 18 populations, and contemporary gene flow by the estimation of recent migration rates among populations. To assess bee-mediated pollen dispersal on a smaller scale, we conducted paternity analyses within a focal population, and quantified within-population spatial genetic structure in four populations. Gene flow was limited to certain nearby populations within continuous forest blocks, whereas drift appeared to dominate on larger scales. Limited long-distance gene flow was supported by within-population patterns; gene flow was biased towards nearby plants, and significant small-scale spatial genetic structure was detected within populations. These findings suggest that, although female euglossine bees might be effective at moving pollen within populations, and perhaps within forest blocks, their contribution to gene flow on the regional scale seems too limited to counteract genetic drift in patchily distributed tropical plants. Among-population gene flow might have been reduced following habitat fragmentation.

Evolutionary consequences of ecological factors: pollinator reliability predicts mating-system traits of a perennial plant.

The reproductive-assurance hypothesis predicts that mating-system traits will evolve towards increased autonomous self-pollination in plant populations experiencing unreliable pollinator service. We tested this long-standing hypothesis by assessing geographic covariation among pollinator reliability, outcrossing rates, heterozygosity and relevant floral traits across populations of Dalechampia scandens in Costa Rica. Mean outcrossing rates ranged from 0.16 to 0.49 across four populations, and covaried with the average rates of pollen arrival on stigmas, a measure of pollinator reliability. Across populations, genetically based differences in herkogamy (anther-stigma distance) were associated with variation in stigmatic pollen loads, outcrossing rates and heterozygosity. These observations are consistent with the hypothesis that, when pollinators are unreliable, floral traits promoting autonomous selfing evolve as a mechanism of reproductive assurance. Extensive covariation between floral traits and mating system among closely related populations further suggests that floral traits influencing mating systems track variation in adaptive optima generated by variation in pollinator reliability.

Contrasting patterns of genetic structuring in natural populations of Arabidopsis lyrata Subsp. petraea across different regions in northern Europe.

Level and partitioning of genetic diversity is expected to vary between contrasting habitats, reflecting differences in strength of ecological and evolutionary processes. Therefore, it is necessary to consider processes acting on different time scales when trying to explain diversity patterns in different parts of species' distributions. To explore how historical and contemporary factors jointly may influence patterns of genetic diversity and population differentiation, we compared genetic composition in the perennial herb Arabidopsis lyrata ssp. petraea from the northernmost parts of its distribution range on Iceland to that previously documented in Scandinavia. Leaf tissue and soil were sampled from ten Icelandic populations of A. lyrata. Seedlings were grown from soil samples, and tissue from above-ground and seed bank individuals were genotyped with 21 microsatellite markers. Seed bank density in Icelandic populations was low but not significantly different from that observed in Norwegian populations. While within-population genetic diversity was relatively high on Iceland (H(E) = 0.35), among-population differentiation was low (F(ST) = 0.10) compared to Norwegian and Swedish populations. Population differentiation was positively associated with geographical distance in both Iceland and Scandinavia, but the strength of this relationship varied between regions. Although topography and a larger distribution range may explain the higher differentiation between mountainous Norwegian relative to lowland populations in Sweden, these factors cannot explain the lower differentiation in Icelandic compared to Swedish populations. We propose that low genetic differentiation among Icelandic populations is not caused by differences in connectivity, but is rather due to large historical effective population sizes. Thus, rather than contemporary processes, historical factors such as survival of Icelandic lineages in northern refugia during the last glacial period may have contributed to the observed pattern.

Genetic constraints predict evolutionary divergence in Dalechampia blossoms.

If genetic constraints are important, then rates and direction of evolution should be related to trait evolvability. Here we use recently developed measures of evolvability to test the genetic constraint hypothesis with quantitative genetic data on floral morphology from the Neotropical vine Dalechampia scandens (Euphorbiaceae). These measures were compared against rates of evolution and patterns of divergence among 24 populations in two species in the D. scandens species complex. We found clear evidence for genetic constraints, particularly among traits that were tightly phenotypically integrated. This relationship between evolvability and evolutionary divergence is puzzling, because the estimated evolvabilities seem too large to constitute real constraints. We suggest that this paradox can be explained by a combination of weak stabilizing selection around moving adaptive optima and small realized evolvabilities relative to the observed additive genetic variance.

Investigation of the genetic diversity of domestic Capra hircus breeds reared within an early goat domestication area in Iran.

BACKGROUND: Iran is an area of particular interest for investigating goat diversity. Archaeological remains indicate early goat domestication (about 10,000 years ago) in the Iranian Zagros Mountains as well as in the high Euphrates valley and southeastern Anatolia. In addition, mitochondrial DNA data of domestic goats and wild ancestors (C. aegagrusor bezoar) suggest a pre-domestication management of wild populations in southern Zagros and central Iranian Plateau. In this study genetic diversity was assessed in seven Iranian native goat breeds, namely Markhoz, Najdi, Taleshi, Khalkhali, Naini, native Abadeh and Turki-Ghashghaei. A total of 317 animals were characterized using 14 microsatellite loci. Two Pakistani goat populations, Pahari and Teddy, were genotyped for comparison. RESULTS: Iranian goats possess a remarkable genetic diversity (average expected heterozygosity of 0.671 across loci, 10.7 alleles per locus) mainly accounted for by the within-breed component (GST = 5.9%). Positive and highly significant FIS values in the Naini, Turki-Ghashghaei, Abadeh and Markhoz breeds indicate some level of inbreeding in these populations. Multivariate analyses cluster Iranian goats into northern, central and western groups, with the western breeds relatively distinct from the others. Pakistani breeds show some relationship with Iranian populations, even if their position is not consistent across analyses. Gene flow was higher within regions (west, north, central) compared to between regions but particularly low between the western and the other two regions, probably due to the isolating topography of the Zagros mountain range. The Turki-Ghashghaei, Najdi and Abadeh breeds are reared in geographic areas where mtDNA provided evidence of early domestication. These breeds are highly variable, located on basal short branches in the neighbor-joining tree, close to the origin of the principal component analysis plot and, although highly admixed, they are quite distinct from those reared on the western side of the Zagros mountain range. CONCLUSIONS: These observations call for further investigation of the nuclear DNA diversity of these breeds within a much wider geographic context to confirm or re-discuss the current hypothesis (based on maternal lineage data) of an almost exclusive contribution of the eastern Anatolian bezoar to the domestic goat gene pool.

Seed dispersal in time can counteract the effect of gene flow between natural populations of Arabidopsis thaliana.

Plants may escape unfavorable environments by dispersing to new sites, or by remaining in an ungerminated state at a given site until environmental conditions become favorable. There is limited evidence regarding the occurrence, interplay and relative importance of dispersal processes in time and space in plant populations. Thirty-six natural populations of the annual ruderal species Arabidopsis thaliana were monitored over five consecutive years, sampling both seed bank and above-ground cohorts. We show that immigration rates are considerably higher than previously inferred, averaging 1.7% per population yr(-1). On the other hand, almost one-third of the individuals in a given above-ground cohort result from seeds shed 2 or 3 yr back in time in 10 of the studied populations. Populations that disappeared one year were recolonized by regeneration from the seed bank the subsequent year. Thus, dispersal in both time and space is an important contributor to the structuring of genetic variability in natural populations of A. thaliana, where a high dispersal rate in time may partly counteract the homogenizing effects of spatial seed and pollen dispersal.

Development of microsatellite markers for the neotropical vine Dalechampia scandens (Euphorbiaceae).

PREMISE OF THE STUDY: Microsatellite markers were developed to assess polymorphism and level of genetic diversity in four Mexican populations of the neotropical vine Dalechampia scandens (Euphorbiaceae). • METHODS AND RESULTS: Thirty-seven microsatellite markers representing bi-, tri-, tetra-, and pentanucleotide microsatellite repeats were developed. In total, 166 alleles were identified across 54 individuals. The number of alleles varied from one to 11 with an average of 4.49 alleles per locus. All loci except one were highly polymorphic between populations, whereas considerably less variation was detected within populations for most loci. The average observed and expected heterozygosities across study populations ranged from 0 to 0.63 and 0 to 0.59, respectively, for individual loci, and a deviation from Hardy-Weinberg equilibrium was observed for most loci. • CONCLUSIONS: The developed markers may be useful for studying genetic structure, parentage analysis, mapping, phylogeography, and cross-amplification in other closely related species of Dalechampia.

Genetic consequences of seed banks in the perennial herb Arabidopsis lyrata subsp. petraea (Brassicaceae).

PREMISE OF THE STUDY: Seed banks may increase the effective population size (N(e)) of plants as a result of elevated coalescence times for alleles residing in the populations. This has been empirically demonstrated in populations of the annual Arabidopsis thaliana, whereas comparable data for perennial species are currently lacking. We studied the contribution of seed banks to effective sizes of natural populations of the self-incompatible, perennial Arabidopsis lyrata subsp. petraea, a close relative of A. thaliana. • METHODS: Fourteen populations of A. lyrata collected throughout the Norwegian distribution range were analyzed using microsatellite markers. • KEY RESULTS: The genetic composition of seed-bank and aboveground cohorts was found to be highly similar, with little genetic differentiation between cohorts in most populations. However, the proportion of private alleles was higher in aboveground than in seed-bank cohorts. The presence of seed banks significantly increased total N(e), but the contribution from seed banks to overall N(e) were lower than the contribution from aboveground cohorts in most populations. Estimated historical N(e) values, reflecting the effective sizes of populations throughout the history of the species, were considerably higher than estimates of contemporary N(e), reflecting number of reproducing individuals within the past few generations. • CONCLUSIONS: The results show that the seed bank contributes to total N(e) in the perennial herb A. lyrata. However, the contribution is similar to or lower than that of the above-ground fraction of the population and markedly weaker than that previously documented in the annual A. thaliana.

Seed banks cause elevated generation times and effective population sizes of Arabidopsis thaliana in northern Europe.

It is commonly found that effective population sizes of natural populations are much smaller than census sizes of plants and animals. However, theoretical studies have shown that factors rarely investigated empirically, like seed banks in plants and diapause in animals, may have profound influence on effective sizes. Here we investigate whether the presence of seed banks can explain the relatively high genetic variability observed in northern European Arabidopsis thaliana populations with small census sizes. We have genotyped three above- and below- ground cohorts in 27 Norwegian populations using single nucleotide polymorphism markers. Although the populations varied extensively in levels of variability within and between cohorts, standard genetic population measures were comparable to those obtained in previous studies on above-ground cohorts using microsatellite markers. Estimated effective population sizes are larger for total populations (containing both seed bank and above-ground cohorts for 1 year) compared to each of the cohorts considered separately. Using a conservative approach, we find that the effective sizes are larger than census sizes of local populations, and that the effective generation time is higher than 1 year (3-4 years, on average), making A. thaliana a perennial semelparous plant at many northern European localities.